vp=1578;
f=@(v)4./pi.^(1/2).*(v.^2./vp.^3).*exp(-(v.^2/vp.^2));
disp('0到vp之间分子的占比')
quad(f,0,vp)
disp('0到3.3vp之间分子的占比')
quad(f,0,3.3*vp)
disp('3e4到3e8之间分子的占比')
quad(f,3e4,3e8)